Real-Time Correlators and Non-Relativistic Holography 



Robert G. Leigh and Nam Nguyen Hoang 

Department of Physics, University of Illinois, 1110 W. Green St., Urbana IL 61801-3080, U.S.A. 



Abstract 

We consider Lorentzian correlation functions in theories with non-relativistic Schrodinger 
symmetry. We employ the method developed by Skenderis and van Rees in which the contour 
in complex time defining a given correlation function is associated holographically with the 
gluing together of Euclidean and Lorentzian patches of spacetimes. This formalism extends 
appropriately to geometries with Schrodinger isometry. 
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1 Introduction 



Correlation functions of operators in strongly coupled conformal field theories can often be com- 
puted using the AdS/CFT correspondence. Euclidean correlators have a long history [U [2] while 
the rich analytic structure of various Lorentzian signature correlators can also be obtained. The 
earliest proposal for the latter was by Son and Starinets[3], and there have also been several elabo- 
rations of that method (see for example [HIS]). Recently, Skenderis and van Rees[6l[7] showed how 
the complex time contour of an arbitrary correlation function can systematically be accounted for 
by gluing together manifolds of various signatures, carefully matching fields at the interfaces. This 
method was used to calculate scalar two-point functions in AdS space, and in asymptotically AdS 
spaces. 

The extension of gauge-gravity duality ideas to spacetimes of Galilean isometrics and field theories 
with non-relativistic invariance [Oj [10] has been of much interest in the recent literature. In par- 
ticular, it is expected that such systems are of more direct relevance to condensed matter models. 
Correlation functions have recently been computed using standard holographic methods for scalars 
P[ini[I5 and for fermions [12] . 

In this paper, we reconsider Lorentzian correlators of non-relativistic systems by directly calculating 
them using the techniques of Refs. [6l [7] in Schrodinger geometries. We consider the time-ordered 
correlator and the Wightman function, as well as thermal correlators. 

2 The Schrodinger Geometry and Scalar Fields 

We consider the d + 3 dimensional Lorentzian geometry [9] 



where z > and 6, L are length scales. This geometry has Schrodinger isometry with dynamical 
exponent equal to two. The Killing vectors are of the form 
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Figure 1: Contours corresponding to the time-ordered correlator and the Wightman function, respectively. 
N is central, and D, H, C form an SL(2, R) algebra. 

Consider a massive complex scalar propagating on the non-relativistic (Lorentzian) geometry with 
action 



{t, x) coordinates-^ is not geometric in the usual sense. The isometry N : ^ + a is central and 
thus N is strictly conserved. Each operator of the boundary theory can be taken to have a fixed 
momentum ('particle number') conjugate to ^. ^ is usually taken compact (with circumference R) 
so that the spectrum of possible momenta is discrete. In this case, the dimensionless ratio b/R is 
a parameter of the theory. 

For example, the graviton mode coupling to the stress energy tensor of the boundary theory has 
particle number zero \13\ I14j. Here, we will consider a complex scalar with definite but arbitrary 
particle number n. As we will see, it is very important that the scalar be complex. First, it carries 
a charge under and so we should expect it to be complex. More importantly though, it is dual 
to an operator in a non-relativistic theory, and in such a theory there is a sort of polarization: a 
simple example of this occurs in free field theories, in which the elementary field creates a particle 
(and not anti-particle) state. 

Now, in this paper we consider correlators of various types. In this regard, as developed by Skenderis 
and van Rees|6|[7|, we regard the metric ([1]) as defined formally for complex t, and a given correlator 
is constructed from a particular contour in the complex t plane. Here, we consider two such cases, 
in which the contour is constructed from horizontal (Lorentzian time) and vertical (Euclidean time) 
contour segments (see Fig. [1]). In the next two subsections, we consider scalar fields in Lorentzian 
time and in Euclidean time, respectively. 




(8) 



The usual interpretation is that the dual theory lives on M^''^ at 2: = and is coordinatised by the 
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2.1 Lorentzian signature 



Given the metric ([T|) for real time, the scalar equation of motion takes the form 

z^dl<i) -{d+ l)zd,(p + z^{2dtd^ + df(p) + b^d^cp - mlcj) = 0. (9) 
We look for solutions of the form 

in which case / satisfies 

z^d^J -{d+ l)zdj + z\2ujn - P)f - m^f = 0, (11) 

where = uiq + n^6^. The general solution of (jlip can be written in terms of modified Bessel 
functions as 

/n,^i(^) = A{uJ,k)zl+'K,{qz)+B{uJ,k)z-^+'Uqz) (12) 



with = Y (f + 1)^ + and (7 = = — 2a;n. i^j^ and /i, correspond to non-normalizable 
and normalizable modes, respectively. Their asymptotic behavior is as follows 

zl+^K,{qz^O) = r(z.)|i^ + ... (13) 

Z2+^Uqz^0) = — + ... (14) 



zt+ii^,(|gz| ^ 00) = J^—e-'^' + ... (15) 



z^+H,{\qz\ ^ 00) = Jf— e'^'il + ...)+ 6-'^'-'^^"+' + ... (16) 

y 27rg L J 

For (7^ < 0, both K,^ and Ii, are regular everywhere, while for q^ > 0, /i, diverges for large z 
and should be discarded. This situation is very similar to that of a scalar field propagating on 
AdSd+3, where the solution can also be written in terms of modified Bessel functions. In fact this 
similarity is very useful and was employed in Ref. [11] to compute the non-relativistic bulk-to- 
boundary propagator. We note though that there is a small but important difference due to the 
non-relativistic nature of the boundary theory, that we will explain presently. 



Without loss of generality, we take n > 0. To construct the most general solution (with fixed n), 
we must integrate over all values of uj,k. However, q has a branch point at to = k'^/2n, and we 
must then say how to integrate over u. Following [6], we do so by moving the branch point off of 
the real lv axis by defining q^ = \/—2ivn + k"^ — ie, q^ = \/—2ivn + fc^ + ie. The branch cut is taken 
along the negative real axis. Clearly, we have made a choice here, but we will see later that this 
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is the correct choice, for physical reasons. Notice that since Re{qe) ■, Re{qe) > 0, Ki, always decays 
exponentially as \qz\ oo. In contrast, the large z behavior of tells us that q,q cannot have a 
real part. As a result, the ie insertion should not be applied for the normalizable mode.^ 



With these comments, we arrive at the general solution to (jlip in Lorentzian signature 
where we have used Iu{y/q^z) = Iu{—i\q\z) ~ Ju{\q\z). 



2.2 Euclidean signature 

Next, we consider a similar analysis in Euclidean signature. To do so, we Wick rotate the metric 
© to[l5] 

ds^ = (b^^ + -^^d-dC+f+dz-^ (18) 

Although this metric is complex and thus not physical, it is possible to trace carefully through the 
analysis, and this is what we need to do in any case for Euclidean signature. 

The general solution is 



(n) , 



,^(„)(r,f) = e— «y'^^e----Hf+il(c^i,,fe)i^.(g-^z) (20) 



where now qE = y Qe — y ^'^ ~ i'^'^En. Note that in this case, the branch point is at imaginary 
UJE-, and so no ie insertion is necessary. 

In contrast to the Lorentzian case, the Euclidean scalar does not have a normalizable mode. This 
is because qE and qE cannot be pure imaginary, so I^iqEz) is never regular in the interior. It is 
important to note, however, that this statement applies to the case r € (— cxd, oo). If r is restricted, 
a normalizable mode can emerge. For example, if r € [0, cxd), we write uje = —i^ for (f) and oje = 
for (j) and the following mode is allowable 

<P ~ e^'^^e-'^^+^^-^zi+^I^(gz) (21) 
~ e-™«e-'^"-^^-^zi+i/^(gz) (22) 

as long as ciJ > and —2ujn + A;^ < 0, or equivalently uj > J? /2n. 



^This fact was not clearly spelled out in Ref. 6 in the relativistic analogue, but we will see later that it is an 
important point. 
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A similar result pertains in the finite temperature case where r S [0, Observe however that in 
contrast to the relativistic real-time formalism, there is no normalizable mode for the Euclidean 
segment if we restrict r € (— cxd, 0). This is because we would need both cj < and —2ujn + A;^ < 0, 
and these contradict each other. This will have important consequences. In particular we note that 
there is no normalizable mode in the segment Mq of either contour in Fig. [TJ 



3 Non-Relativistic Holography and Correlators 
3.1 Matching Conditions 

To construct correlation functions, we must match solutions at the interfaces between contour 
segments. We will label field values on a contour segment M„ by a subscript, (/>„. Let us begin 
by considering the Lorentzian(Mi)-Lorentzian(M2) interface in Fig. [1)3, where ti G [0,T] and 
t2 G [T, 2T] (where T — > 00 is a large time). The total action (for these two segments) is 

T , _ _ f-2T 

dti [g^j^^^d^(t)idy(l)i + ml /L^ 4)1(1)1 j - J 



S = Smi + = I dti {g^f^i^d^(l)idy4)i + ml/L'^(t>i4>i^ - j dt2 (^gM^d^4>2d,y(l)2 + ml / L'^ (t)2ct>2^ 

(23) 



The relative minus sign arises because Mi and M2 have opposite orientation. For the same reason, 
the metric in M2 is 

dsl,^ =L^{-§+ -^dt2dC+f + dz^ ^^ , (24) 
which has an extra minus sign in the off-diagonal component. 

Requiring continuity of the momentum conjugate to (j) at the intersection ti = ^2 = T, we get 

d^<pi = d^cp2. (25) 
Along with the continuity of we conclude that the matching conditions at ti = t2 = T are 

cl^iiT) = <j)2{T) (26) 

Til = 71-2 (27) 

Thus, we do not need to impose first-order time derivative continuity of fields along the contour 
as in the relativistic case — it is just replaced by particle number conservation. It turns out that 
()26|27p are also the matching conditions for Euclidean - Lorentzian interfaces. 
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3.2 Convergence and the Choice of Vacuum 



The non-relativistic holographic correspondence is in general the same as its relativistic counterpart, 
where the path integral with specified boundary conditions in the bulk is identified with the partition 
function with sources inserted in the boundary theory. In the case of a complex bulk scalar, we 
must temporarily treat the sources 0(o) and 0(o) independent. The near boundary expansion of 
the fields are qualitatively the same as scalars on AdSd+3 

<Pin) = e'<{{z^- (,^(0) + zV{2) + o{z')) + (t;(o) + z\,) + o{z^)) } (28) 
^(n) = e'^^[{z^- {4>^o) + ^'<^(2) + o{z^)) + z^+ (i;(o) + z\,) + o{z^)) }, (29) 
with A± = 1 + d/2 lb 1^ and 

"^(2-) = 2m(2A+-(d + 2)-2m)°°'^(2m-2) , (30) 

where here Dq = 2indt + is the non-relativistic Laplacian. As usual the holographic correspon- 
dence implies 



JS', 



bulk r 



where C denotes the contour. Although we have a very different geometry, it's easily seen that in 
each patch of the contour the bulk (either Euclidean or Lorentzian) on-shell action 

Sos = \ ^d''^^xdi^\ ^ g'' d,4> (32) 

is essentially the same as scalars on AdS^+z- As a result, the renormalization procedure proceeds 
in the same way as AdSd+'i/CFT(i+2-, which was carried out in much details in [8]. In specific, for 
Lorentzian signature the counter terms take the form, 

Set = jd''+^xdi^[ ^^'^~^^ 4><l^ + +...), (33) 

where the {d + 2)-dimensional induced metric determinant and D-y = z^{2indt + 

df) (we will set L = 1 from now on). The dots represent higher derivative terms. For special cases 
where v is an integer, logarithmic counter terms ~ loge may appear It's important to note 
that Set preserves the Galilean subalgebra, since |0^,i^j] = 0. This is in parallel with relativistic 
holography where the Poincare subalgebra is preserved by the counter terms. In any case, v^q) will 
determine the v.e.v of the dual operator and its derivative with respect to the source 0(o) gives us 
the 2-point functions. 

There is, however, a subtlety of which we must be cognizant. Unlike relativistic field theories, in 
non-relativistic field theories an elementary field ^ and its Hermitian conjugate play the role 
of creation and annihilation operators. There is a freedom to choose which is an annihilator, or 
equivalently a freedom to pick the vacuum. Once a convention is chosen, ^ and ^'^ are no longer 
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on the same footing. This is also true for any operator O, 0\ in which O is constructed only from 
annihilators. This corresponds to the fact that there is only a single pole in the complex w-plane in 
the non-relativistic case. Consequently, the time-ordered propagator will in fact have only a single 
temporal ^-function present. We expect to see this coming about in the analysis, but to see this 
properly, one has to be careful with the convergence of various integrals. 



4 Correlation Functions 



In both cases shown in Fig. [H we have an initial vertical contour Mq. The correlation functions 
of interest are computed by including source(s) on horizontal component (s) of the contour. We 
first show that given such a contour component Mq, there is no normalizable mode (such a mode 
would be everywhere subleading in the z — > expansion). This implies that any solution with a 
specific boundary condition is unique. Indeed, we argued in Section [2^2] that there is no non-trivial 
normalizable solution in Mq. So in the cases of interest (no sources on Mq), i;^o = identically. The 
matching condition between (pQ and (pi then requires that 01(^1 = Q,x,z) = 0. The most general 
normalizable solution on Mi is 

rr'^{h,x,z) = j e— *^+^Hf+ie(-g2)i?(^,fc)j,(|g|^). (34) 

Multiply by e~'^'^^~'^^' Jy{\q'\z) with q'"^ = —2Lo'n + k''^ < and integrate over x and z. We 
then find 

= I ^^d'^ e^^<^-^^'^B{u,,m-q^)(^j^ dz z.h{\q\z)U\q'\z)) (35) 
The z- integral is elementary (see Appendix, eq. (I59p ) and this becomes 

= y'^^d'^xe-(^^-^^')i?(..,fc)0(-52)^5(|^|_|^/|) (3g) 

= - [ —B{uj, k')e{2ujn - k'^)S{uj - Lo') (37) 
n J 2tt 

= ^B{J,k')e{-q'^). (38) 

Thus, if (j)i{t,x,z) = at some time, there is no non-trivial normalizable mode. This reasoning in 
fact applies for all segments of both contours in Fig. [H 



4.1 Bulk-Boundary Propagator and Time-ordered Correlator 

Given the absence of a normalizable mode, any solution with sources that we find for the two 
contours in Fig. [T] is unique. In this subsection, we consider contour Fig. [T^, with segments 
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Figure 2: Contour of integration in the complex p-plane for the Lorentzian bulk-boundary propagator. 



(to G (-00,0]), Ml (ti G [0,r]), M2 {t2 G [0,00)). We place a single (5-function source at 
X = Q^ti = ti on Ml. From our discussions above, 0i must be of the form 



(t>l,{n){tl,X,z) 



Ky{q^z). 



(39) 



as this satisfies 2;~^-</)i („)(ti, x, z) [^^^^ = e*"^5(ti — ti)(5(x), and any ambiguity corresponds to 
normalizable modes, which we have argued are zero. Since there are no sources on M2, (j)2 takes 
the form 



52, (n) 



2-K (27r)' 



-u}(T+iT—iti )+ik-x 



Jui\Q\z)- 



(40) 



which has been deduced from the matching condition (/>i(ti = T) = 4'2{t = 0) as follows. For any 
time ti > ti, we can re-expand (pi in terms of Ju's. In particular, at ti = T, we should have 



duj d'^k 
27r (27rp 



duj d'^k 
2^(2^ 



-i-^{T-t,)+-ik-Xfj^^^ k)e{-q^)zJMz) (41) 



for some C(a;, k). To find this coefficient we use the same trick as in the last subsection: multiply 
both sides by e^'^' Ju{\q'\z) with q'"^ = —2uj'n + k'"^ < and integrate over x,z. The 
right-hand side gives Y-^e{-q'^)C{uj',k'), while the left-hand side can be computed using (j60p to 



give ^I^T- 



The bulk-boundary propagator is essentially identified with (f>i itself: if we simply strip off the e™^ 
factor, we can write 



Kn,n'{t,X,z) = (5„,„/K(„)(t,X,z) 

, ^ . 2zi+'^/2 f duj d'^k 

K(n){i,x,z;t) = 



^iu;(t^i)+^k■x f Qe\ ^^^^^ 



T{iy) J 27r {2ttY ' V 2 



(42) 
(43) 



As shown in Ref. [11] for example, this is closely related to the bulk-boundary propagator in 
AdSd+3- Alternatively, we may perform the integration directly, following the analogous treatment 
in Ref. [6]. To do so, it is convenient to convert the w-integral to an integration over p = qe, and 
the contour in the p-plane is as shown in Fig. [21 Here though there is just one branch point (at 
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LO = k'^/2n — ie) and the ie tells us in which sense to traverse the cut. One arrives at 

A-,„,(M-..;*) = «fe-t,)^^y [j^J e (44) 

where A± = 1 + d/2 ± u. 

The correlator is then identified with the z'^+ coefficient in the near boundary expansion of 4>i 
(without the e*"^ factor) 

(r(6,,(.,.)o;„,(^^.*'.))) - ^ iT" (^'•"^' (-) 

4.2 Wightman function 

The time-ordered correlator, as we have explained, contains a single temporal ^-function. It does 
not tell us about {0{x,ti)0^ [x' ,t[)) for t[ > ti. To find this 2-point function we work with the 
contour of Fig. Denote the segments by Mq (tq G (-oo,0]), Mi {h E [0,r]), M2 (t2 e [r,2r]) 
and M3 (t3 E [0, 00)) as sketched in the figure. We place a 5-function source at x = 0,ti = ti on 
Afi and nowhere else. The Wightman function is obtained then from <j)2, the field on M2. Here 
(pQ = and (pi remain the same as ()39p . Given experience from the last subsection, we can see 
immediately that (/>2 should be 

^..H = ^/^'^'^'^'l e-(— *^)+--(M)%(_,^)j,(|,|,). (46) 

This has been determined by requiring the matching condition (/>i(ti = T) = (j)2{t2 = T). Notice 
the unusual e+*'^*2+«fc a; ^^ve factor. It is related to the fact mentioned before that along this part 
of the contour, the metric has an extra minus sign in the gt^^^ component. 

It is now necessary to compute (j)2 in coordinate space. We make a change of variable p = \q\ = 



^ nT{v)2'' Jo ^ -v^ ^27r) 

We note that both integrals converge if 2T — t2 — ti ^ 2T — t2 — ti — ie. The first integral can be 
computed using ([6T]) . while the second one is just a Gaussian integral. The final result is 



1 /n\A+-i / z \A+ -2+5" 



where ^2 = 2T — t2- Observe that (/)2 is closely related to the bulk-boundary propagator except 
for the absence of the step function and a different ie insertion, as expected. 
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Figure 3: Thermal contour. Points with a circle are identified. 
The vacuum expectation value of (D{i2,x) is 



m-^ — 



Taking a derivative with respect to (pK^o) and setting the source to zero, we get the Wightman 
function 



m—^ — 



(0(t2,x)0t(ti,x')) = ^^7^ (^) ^ (50) 

7r"/"T(z^) (^2 — ii — ^e) + 

Notice that is always in the front of O because ti is always the earlier contour time. 



4.3 Thermal Correlator 

Finally, we compute a thermal correlator by taking the time direction to be compact of period (3. 

To compute the thermal time-ordered correlator and Wightman function, we consider the thermal 
contour shown in Fig. [3l where t = and t = — i/3 are identified. We place a (5-function source at 
ti = = 0. Note that in contrast to the previous discussions, here there is no Mq component of 
the contour. It is convenient in this context to write the general solution along Mi in the form 



^ (51) 
where q-e = Qe = \l —lion + fc^ + ie. In order that this correspond to a 5-function source for 
z — > 0, we must have B = \. (Furthermore, the case B = —A corresponds to a normalizable 
mode.) Note that because of the condition on A, B, although A and B are not necessarily analytic 
functions, their sum is analytic. Thus for example, for any pole in A, there will be a corresponding 
pole in B with opposite residue. All of their poles will contribute opposite residues and cancel out 
each other in the limit e — > 0. In (j5ip . the first term has support for ti > ti, while the second has 
support for ti < ti. 
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The matching condition at (Mi,M2) and (M2,M3) intersections imply that 



T{u) J 27r {2ttY 



'iuj{2T-t2-ti)+ik-x 



-uj{T-j—iti)+ik-x 



Aiu;,k)(^] M\q\z)ei-q') (52) 



(53) 



The thermal condition 0i(ti = 0) = </>3(t3 = P) along with A + B = 1 then gives 

1 1 



A 



B 



As usual, the time-ordered propagator is the coefficient of in the small z expansion of (j)i 
(without the e*"^ factor). Hence we get^ 



{T(d{x)d\x') 



duj d'^k 



^~Mt-t')+ik-{x-x') ^ {-2ujn + k'^ -ief ^ {-2u;n + P + i^y 



1 - 6-/3'-' 



2tt {2-kY 

Note that this has the expected form for a thermal correlator^ 

(r(o(x)Ot(x'))) = -iV(c^)AA(c^,fc) + {l + N{u:))^R{u:,k) 



1 - eP"^ 



(55) 



(56) 

In the present notation, N = —B. We can also write this as the zero temperature result plus a 
finite temperature piece: 



(rfo(x)Ot(x') 



duj d'^k 
2tt {2'kY 



-iui{t-t')+ik-{x-x') 



J2v 



1 



1 - eP' 



2u „2us 



(57) 



The Wightman function can also be read off from (j)2 

duj d'^k 



(O(x)O^(x')) 



27r {2T:y 



^ i^it-tf-r^+i±ix-x') (2^^ - ^y e{2u:n - P) 



1 - e-^"^ 



(58) 



A Appendix 



We record integrals that have been useful in the above analysis. 



f-OO 2 

/ tJ^{qt)Mq't)dt=-5{q-q') 

Jo q 



q,q' veal, > 



K^{at)Mbt)t''+''+^dt 



i2aY{2bYT{fi + iy+l) 
(a2 + 62)M+^+l 

Re{iy + l) > Re{fi), Re{a) > \Im{b)\ 



f 

Jo 



e-" * t''+^Mbt)dt 



b"" 



{2a 



i?e(i/) > -1, i?e(a2) > 



^For integer v, there is an extra logarithmic factor, namely is replaced by q±^ lng±£. 



(59) 



(60) 
(61) 
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